On the electric conductivity of hot QCD matter 
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We study the electric conductivity of hot QCD matter at various temperatures T within the off- 
shell Parton-Hadron-String Dynamics (PHSD) transport approach for interacting partonic/hadronic 
systems in a finite box with periodic boundary conditions. The response of the strongly-interacting 
system in equilibrium to an external electric field defines the electric conductivity ao- We find a 
sizeable temperature dependence of the ratio ao/T well in line with calculations in a relaxation time 
approach for Tc<T < 2.5Tc. The ratio drops in the hadronic phase with T", shows a minimum close 
to Tc and becomes approximately constant (^0.3) above ~5Tc. Our findings imply that the QCD 
matter even at T^Tc is a much better electric conductor than Cu or Ag (at room temperature). 
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High energy heavy-ion reactions are studied experi- 
mentally and theoretically to obtain information about 
the properties of nuclear matter under the extreme condi- 
tions of high baryon density and/or temperature. Ultra- 
relativistic heavy-ion collisions at the Relativistic Heavy- 
Ion Collider (RHIC) and the Large Hadron Collider 
(LHC) at CERN have produced a new state of matter, 
the strongly interacting quark-gluon plasma (sQGP), for 
a couple of fm/c in volumes up to a few thousand frrt' in 
central reactions. The produced QGP shows features of a 
strongly-interacting fluid unlike a weakly-interacting par- 
ton gas [l| as had been expected from perturbative QCD 
(pQCD). Large values of the observed azimuthal asym- 
metry of charged particles in momentum space [5] could 
quantitatively be well described by ideal hydrodynam- 
ics up to transverse momenta of 1.5 GeV/c [3|. Recent 
studies of 'QCD matter' in equilibrium - using lattice 
QCD calculations [4, 5] or partonic transport models in 
a finite box with periodic boundary conditions |.6j, [Tj - 
have demonstrated that the ratio of the shear viscosity 
to entropy density rj/s should have a minimum close to 
the critical temperature Tc, similar to atomic and molec- 
ular systems [J |^ . On the other hand, the ratio of the 
bulk viscosity to the entropy density C,/s should have a 
maximum close to Tc 0] or might even diverge at Tc [lol | . 
Indeed, the minimum of r;/s at « 160 MeV is close to 
the lower bound of a perfect fluid with rj/s = l/(47r) [llj 
for infinitely coupled supersymmetric Yang-Mills gauge 
theory (based on the AdS/CFT duality conjecture). This 
suggests the 'hot QCD matter' to be the 'most perfect 
fluid' [1, [H, [Ts'l . On the empirical side, relativistic vis- 
cous hydrodynamic calculations (using the Israel-Stewart 
framework) also require a very small rj/s of 0.08 — 0.24 



in order to reproduce the RHIC elliptic flow V2 data [IJ] ; 
these phenomenological findings thus are in accord with 
the theoretical studies for ?]/s in [tI [Tsl [l6j. 
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Whereas shear and bulk viscosities of hot QCD matter 
at finite temperature T presently are roughly known, the 
electric conductivity ao is a further macroscopic quan- 
tity of interest |17j|. The basic question is: Is the 'hot 
QCD matter' a good electric conductor? At first glance 
one might expect the deconfined QCD medium to be 
highly conductive, since color charges - and associated 
electric charges of the fermions - might move rather 
freely in the colored plasma. However, due to the ac- 
tual high interaction rates in the plasma - reflected in 
a low ratio rj/s - this expectation is not so obvious. 
First results from lattice calculation on the electromag- 
netic correlator provide results that vary by more than 
an order of magnitude [T8l - [2^ . Furthermore, the con- 
ductivity dependence on the temperature T (at T>Tc) 
is widely unknown, too. The electric conductivity ctq is 
also important for the creation of electromagnetic fields 
in ultra-relativistic nucleus-nucleus collisions from par- 
tonic degrees-of-freedom, since ao specifies the imaginary 
part of the electromagnetic (retarded) propagator and 
leads to an exponential decay of the propagator in time 
^ exp{—ao{t — t') / {he)) High values of ao would thus 
lead to the screening of external electromagnetic fields 
in the bulk of the highly-conducting quark-gluon plasma 
similar to the Meissner effect in super-conductors as well 
as the "skin-effect" for the electric current. Accordingly, 
a sufficient knowledge of ctq {T) is mandatory to explore a 
possible generation of the Chiral-Magnetic-Effect (CME) 
in predominantly peripheral heavy- ion reactions |24| . 

In this work we extract the electric conductivity 
ao(T) for 'infinite parton/hadron matter' employing the 
Parton-Hadron-String Dynamics (PHSD) transport ap- 
proach [2^ , which is based on generalized transport equa- 
tions derived from the off-shell Kadanoff-Baym equa- 
tions [26j for Green's functions in phase-space represen- 
tation (beyond the quasiparticle approximation). This 
approach describes the full evolution of a relativistic 
heavy-ion collision from the initial hard scatterings and 
string formation through the dynamical deconfinement 
phase transition to the strongly-interacting quark-gluon 



plasma (sQGP) as well as hadronization and the subse- 
quent interactions in the expanding hadronic phase. In 
the hadronic sector PHSD is equivalent to the Hadron- 
String-Dynamics (HSD) transport approach [27*1 ~ a co- 
variant extension of the Boltzmann-Uehling-Uhlenbeck 
(BUU) approach (ill - that has been used for the de- 
scription of pA and AA collisions from lower SIS to RHIC 
energies in the past. On the other hand, the partonic 
dynamics in PHSD is based on the Dynamical Quasi- 
Particle Model (DQPM) [H, which describes QCD prop- 
erties in terms of single-particle Green's functions (in the 
sense of a two-particle irreducible (2 PI) approach) and 
reproduces lattice QCD results - including the partonic 
equation of state - in thermodynamic equilibrium. For 
further details on the PHSD off-shell transport approach 
and hadronization we refer the reader to [a, [2^ [30, HH . 

Here, we concentrate on calculating the electric con- 
ductivity for 'infinite' QCD matter, which we simulate 
within a cubic box with periodic boundary conditions 
at various values for the energy density (or tempera- 
ture). The size of the box is fixed to = 9^ fm"^ as 
in our previous investigations 0, 0] ■ The initialization is 
done by populating the box with light (u, d) and strange 
(s) quarks, antiquarks and gluons slightly out of equi- 
librium. The system approaches kinetic and chemical 
equilibrium during its evolution within PHSD. For more 
details on the simulation of equilibrated partonic systems 
using PHSD in the box we refer the reader to Q. 

Let us remind that PHSD is an off-shell transport ap- 
proach that propagates quasi-particles with broad spec- 
tral functions. Numerically, the continuous spectral dis- 
tribution of mass of a particle (given by its spectral func- 
tion) is probed by a large number of test-particles with 
(evolving) masses Mj(t). In order to include the effects 
from a constant external electric field E^, the propaga- 
tion of each test-particle (in z— direction) is performed 
with the additional force in the equation of motion: 



(1) 



where qj denotes the fractional charges of the test- 
particles (±1/3, ±2/3). The electric current density jz{t) 
then is given by 



eqj 



(2) 



where Mj{t) is the mass of the test-particle j at time t. 
A note of caution has to be given, since due to an ex- 
ternal field we deal with an open system with increasing 
energy density (temperature) in time. Therefore we em- 
ploy sufficiently small external fields, such that the en- 
ergy increase during the computation time stays below 
2% and the increase in temperature below 1 MeV. 

Fig. [1] presents the time dependence of the electric 
current jzit) induced by an external (constant) electro- 
magnetic field of strengths eE^ = 5 MeV/fm, eE^ = 
10 MeV/fm and eE^ = 20 MeV/fm (for the temperature 
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FIG. 1: (color online) The electric current density jz{t) in- 
duced by an external (constant) electric field of strength eE^ 
as a function of time t. The calculations are performed for a 
system of partons at temperature T — 190 MeV in a box of 
volume 729 fm'^ within the PHSD approach. 

T = 190 MeV). It is seen that the current achieves an 
equilibrium value (denoted by jeq) that is proportional 
to the external field. In fact, we obtain the conductivity 
from the ratio of the current density and the electric field 
strength 

a£ _ jeq 

T ~ E,T' 



(3) 



which is shown in Fig. [5] as a function of the external 
field eEz- All results are compatible with a constant ratio 
([3]) as indicated by the straight line. We note that our 
numerical results for cto do not depend on the volume V 
of the box within reasonable variations by a factor of 8. 

We have performed the PHSD studies for par- 
tonic/hadronic systems at various temperatures from 
T = 100 MeV up to T = 350 MeV. The respective resuhs 
for the ratio ctq /T versus the scaled temperature T/Tc are 
displayed in Fig.[3]by the full round symbols. We observe 
a decreasing ratio (Tq/T with T/Tc in the hadronic phase, 
a minimum close to Tc and an approximately linear rise 
with T/Tc above (=158 MeV). Within the error bars 
of our calculations (which in Fig. [3] are indicated by the 
size of the symbols above Tc), the conductivity in the 
partonic phase is described by 



0.01 + 0.16^ 



.T-Tc 



(4) 



T T 

for Tc<T<2.2Tc. The IQCD numbers \lmd^\ are repre- 
sented by symbols with error bars (using Cem = 2e^/3, 
= Aira, a = 1/137). In view of the pQCD prediction of 
a constant asymptotic value for on/T 5.9769/e^ « 65 
in leading order of the coupling [3 HI, HI] , such a linear 
rise of the ratio with temperature might be surprising, 
but it can be understood in simple terms as demonstrated 
below. Only at the highest temperatures studied (~ 5Tc) 
we will see a stabilization of the ratio. 

We recall that the electric conductivity of gases, liq- 
uids and solid states is described in the relaxation time 
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FIG. 2: (color online) The ratio of the conductivity to the 
temperature at T= 190 MeV as a function of the external 
electric field eE^. The statistical accuracy of the calculations 
is reflected in the error bars. The straight line gives the best 
fit to the dimensionless ratio ao/T. 



approach by 



(5) 



where n denotes the density of non-localized charges, r is 
the relaxation time of the charge carriers in the medium 
and m* their effective mass. This expression can be di- 
rectly computed for partonic degrees-of-freedom within 
the DQPM, which was used to match in PHSD the quasi- 
particles properties to lattice QCD results in equilibrium 
for the equation-of-state (EoS) as well as various corre- 
lators 29]. We note that the electromagnetic correlator 
from IQCD calculations [l^ appears to match rather well 
the back-to-back dilepton rate from PHSD at T = lAbTc 
(cf. Fig. 2 in jiSj]), which suggests that the results of our 
calculations for ctq should also be close to the IQCD ex- 
trapolations from 

In the DQPM, the relaxation time for 
quarks/ antiquarks is given by r = l/rg(T) [s^ . 
where Tq{T) is the width of the quasiparticle spectral 
function (cf. [H, H^l)- Furthermore, the spectral distri- 
bution for the mass of the quasiparticle has a finite pole 
mass Mq{T) that is also fixed in the DQPM, as well 
as the density of (m, u, d, J, s, s) quarks/antiquarks as a 
function of temperature (cf. [29, |30||). Thus, we obtain 
for the dimensionless ratio ([3]) the expression [32|] 



^o(r) 



9Mg{T)Tg{T)T' 



(6) 



where nq(T) denotes the total density of quarks and anti- 
quarks and the prefactor 2/9 reflects the flavor averaged 
fractional quark charge squared (X]/'Z/)/3- The result 
for the ratio ([6]) is displayed in Fig. [3] (dash-dot line) 
and does not involve any new parameters. Apparently, 
the PHSD results in equilibrium and the relaxation-time 
estimates match well up to ~2Tc, which demonstrates 
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FIG. 3: (color online) The ratio ao/T as a function of the 
scaled temperature T/Tc {Tc = 158 MeV). The fuU round 
symbols show the PHSD results, the solid blue line is the lin- 
ear fit to the PHSD results (above Tc), while the dash-dotted 
red line gives the corresponding ratio in the relaxation-time 
approach (employing the DQPM parameters). The scattered 
symbols with error bars represent the results from lattice 
QCD calculations: triangles - Refs. [13], diamonds - Ref. [l^. 
squares - Ref. [l^], star - Ref. [2l[|, open circle - Ref. . 
We used for the average charge squared Ceai ~ Sttq/S with 
1/137. Note that the pQCD result at leading order be- 



yond leading log 34] is ao/T ^ 5.97/e^ 



65. 



again that PHSD in equilibrium is a proper transport 
realization of the DQPM Q. 

Our results from the DQPM suggest that above T ~ 
5Tc the dimensionless ratio ([6]) becomes approximately 
constant (« 0.3). This comes about as follows: At high 
temperature T the parton density scales as ~ T^, while 
Mq{T) - T and TqiT) ~ T. Accordingly the ratio §^ 
is approximately constant. Note, however, that energy 
densities corresponding to T > 5Tc are not reached in 
present experiments with heavy- ions at RHIC or LHC! 
On the other hand, (Tq/T rises with decreasing temper- 
ature below Tc (in the dominantly hadronic phase), be- 
cause at lower temperatures the system merges to a mod- 
erately interacting system of pious, which in view of Eq. 
(5) has a larger charge (squared) to mass ratio than in 
the partonic phase as well as a longer relaxation time. 

In summary, we have evaluated the electric conduc- 
tivity cfq{T) of the quark-gluon plasma as well as the 
hadronic phase as a function of temperature T by em- 
ploying the Parton-Hadron-String Dynamics (PHSD) off- 
shell transport model in a finite box for the simulation 
of dynamical partonic/hadronic systems in equilibrium. 
The PHSD approach in the partonic sector is based on 
the lattice QCD equation of state of [35|; accordingly, it 
describes the QGP entropy density s{T), the energy den- 
sity e{T) and the pressure p(r) from IQCD [6, 25, 3^ 
very well. Studies of the QCD matter within PHSD have 
previously given reasonable results also for the shear and 
bulk viscosities rj and C versus T We find in the 
present study that the dimensionless ratio ao/T rises 
above Tc approximately linearly with T up to T = 2.5Tc, 



4 



but approaches a constant above 5Tc, as expected from 
pQCD. This finding is naturaUy explained within the 
relaxation-time approach using the DQPM spectral func- 
tions. Below Tc the ratio ctq/T rises with decreasing tem- 
perature because the system merges to a moderately in- 
teracting gas of pions with a larger charge to mass ratio 
than in the partonic phase and a longer relaxation time. 

The actual values for the electric conductivity o'o(T) 
show that the sQGP even at its minimum (at T « Tc) 
is a better electric conductor than Cu or Ag (at room 



temperature) by about a factor of 500 (using 27rh/e'^ = 
2.58 • lO^ri [3g|). Furthermore, the damping of the elec- 
tromagnetic propagator in the partonic medium (T) — 
{hc)/ao{T) is moderate for temperatures T < 2Tc, yet 
becomes shorter than 2 fm/c for T > 2.8Tc. This sug- 
gests the potential importance of the 'skin-effect' in the 
response of the partonic medium as created in the mid- 
rapidity region of heavy-ion collisions at the LHC to the 
electromagnetic fields generated by the spectators. 
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